Dissipative structures in optomechanical cavities 
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We analyze the possibility of pattern formation in an optical cavity in which one of its mirrors 
can be deformed by radiation pressure. Our model treats the deformable mirror as tense membrane 
that can oscillate, not only back and forth around its axial equilibrium position (center of mass 
motion), but also in its transverse degrees of freedom, vibrating like a drum in combinations of 
modes. We demonstrate the existence of periodic patterns and localized structures (cavity solitons) 
in this model, which should be realizable with current technology. 
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Introduction. Since the eighties of the past century, 
the advent and rapid growing of the field of quantum 
information have boosted the research on quantum tech- 
nologies that try to develop nano and microdevices allow- 
ing for robust implementations of controllable quantum 
interactions, together with long coherence times. By to- 
day the state of the art permits the fabrication of devices 
functioning in such strong coupling regime with unprece- 
dented control and accuracy. These new quantum de- 
vices arc very diverse including cavity-QED [H, optical 
lattices , trapped ions Q , superconducting circuits [if , 
quantum dots [l| , atomic ensembles @ , optomechanical 
devices Q, etc. 

In optomechanical systems Q the interaction occurs 
between a light field and a mechanical oscillator via radi- 
ation pressure. These type of systems have been known 
for a long time in, say, macroscopic optics H, and to 
some extent resemble a nonlinear Kerr cavity |9| as the 
cavity length (and consequently the cavity resonance) de- 
pends on the intracavity field intensity. At present, the 
stress is made on their capability to show quantum co- 
herent phenomena such as cooling and amplification flo| , 
or strong (linear) cou pling effects like optomechanically 
induced transparency [ll|. Attention has also been re- 
cently pai d to the nonlinear dynamics of optomechanical 
arrays [12J , cavities in which radiation pressure competes 
with the photothermal effect [l3| . planar dual-nanoweb 
waveguides subject to radiation pressure [l4[, an d op - 
tomechanical cavities containing atomic ensembles [151 ] . 

Except for some of these recent works, up to now most 
of the studies have considered a small number of modes 
either in the radiation or the mechanical fields. How- 
ever multimodc configurations can provide optomechani- 
cal systems with new capabilities, as the nonlinear inter- 
play between many optical and mechanical modes entails 
the existence of correlations between them. At the quan- 
tum level, these correlations may lead to multipartite 
entanglement At the classical level, they help 

for the spontaneous appearance of dissipative structures 
that are long range ordered field configurations, including 



periodic and quasiperiodic patterns as well as localized 
structures (l9| . which may even exhibit nontrivial tem- 
poral behavior. 

The interplay between the quantum and classical per- 
spectives in multimode systems has received theoreti- 
cal attention in the past mainly in the context of opti- 
cal parametric oscillators, and some exciting phenomena 
have been predicted [2(| . However no such studies exist 
concerning optomechanical systems. 

Here we address the modeling of pattern formation in 
an optomechanical oscillator. We will keep our study at 
the classical level by concentrating on the analysis of the 
general conditions under which nonlinear patterns can 
appear and how these patterns are, and leave the study 
of their quantum properties for a future publication. As 
we show below, periodic patterns and localized structures 
(cavity solitons) are predicted in both 1 and 2 transverse 
dimensions in a wide region of the parameter space. 

The model. Consider an optical cavity with large- 
area mirrors one of which is plane, partially transmitting, 
and immune to radiation pressure because of its stiffness 
and mass, while the other is a perfectly reflecting de- 
formable mirror. The latter can be thought of as a tense 
membrane that can oscillate, forced by the intracavity 
radiation pressure, not only back and forth around its 
axial equilibrium position (center of mass motion) but 
also in its transverse degrees of freedom, vibrating like 
a drum in combinations of modes. The movement of 
the center of mass will be modeled as a homogeneous 
mode (independent of the transverse coordinates) while 
the modes deforming the flatness of the mirror will be 
modeled as the nondispersive surface waves occurring in 
a tense membrane. Energy is fed in the cavity from the 
outside by injecting a coherent laser beam through the 
partially transmitting mirror. We denote by z = 0, L the 
planes containing the coupling mirror and the surface of 
the deformable mirror at rest (i.e., in the absence of il- 
lumination) , respectively. Find a sketch of the system in 
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FIG. 1: Cartoon of the optomechanical system studied in this 
Letter. 



wc choose to be normalized as 



d 2 r u-j (r)u z * (r) = Sjimj/a, 



(4) 



where S is the mirror surface, a is its surface mass den- 
sity, and rrij is the so-called effective mass 23[ which 
won't play any relevant role in our case. As usual [23| . 
the displacement field is expanded onto the normal mode 
basis as 



Q(r,i)=5>;(t)%(r) 



(5) 



mirror will be considered to be a paraxial, coherent beam 
of the form 

E ini (z, r, t) = iVA„j (z, r, t) <?(.*<>-<*>*) + c.c, (1) 

where r = (x, y) denotes the position in the plane trans- 
verse to the cavity axis (z-axis), and V is a constant 
having the dimensions of voltage. The intracavity field 
E (z, r, t) can be written generically as 



E(z,r,t) = iV (A+e 



A-e- lk ^ z ) e - *" L * + c.c. 



(2) 



which is the superposition of two waves of complex am- 
plitudes A± (z, r, t), propagating along the positive (A + ) 
and negative ( A_ ) z direction. Our first goal is to derive 
an evolution equation for the amplitude A + (z = L,r,t) 
at the deformable mirror surface, which we will denote by 
A{r,t) [2l[. The procedure follows the lines of [22| and 
consists in propagating the complex amplitude along a 
full cavity roundtrip, which, under reasonable assump- 
tions (see the supplemental material for details) leads to 



d t A 



-l + iA + iL 2 V' 



where j c is the cavity damping rate, A = (wl — uj c ) /j c 
is the dimcnsionless detuning parameter (uj c is the fre- 
quency of the longitudinal cavity mode closest to wl), 
l c is the diffraction length, = d 2 + d 2 is the trans- 
verse Laplacian, Q (r, t) is the displacement of the mirror 
perpendicular to its flat state (Q = at rest), T is the 
transmissivity of the fixed mirror, and 8 (r, t) is a scaled 
version of the injection field amplitude. 

The equation for the light field must be complemented 
with an equation for the displacement field Q(r, t). In 
the absence of forcing and dissipation Q obeys a wave 
equation d 2 Q + hQ = 0, where L is a suitable differen- 
tial operator. It is useful to diagonalize the problem by 
introducing a set of normal modes {uj (r)}°^_ , which are 

the solutions to the linear problem huj (r) = £l 2 Uj (r) 
subject to appropriate boundary conditions; here Qj is 
the mode j frequency. Without loss of generality we take 
the normal modes as dimcnsionless functions of r, which 



and then the generalized coordinates qj are modeled as 
damped harmonic oscillators forced by radiation pres- 
sure: 



(0) 



In this equation jj is a damping rate and the radiation 
pressure force can be written as 



Fj (*) 



2hk c 
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d 2 r\A(r,t)\ 2 u* (r) 



(7) 



with t c = 2L/c the cavity roundtrip time. In writing Fj 
we used the expression for the radiation pressure on the 
movable mirror P ra( j (r, t) = (5+ (r, t)) /c in terms of the 
time averaged modulus of the Poynting vector associated 
to the wave impinging the deformable mirror, which in 
our notation reads (5+) = 2e cV 2 \A\ 2 [Hj]. 

Up to here the treatment is general. In order to ob- 
tain an evolution equation for Q (r, t) a specific model 
for the the mirror dynamics must be chosen, both for 
its Hamiltonian part, governed by the differential oper- 
ator L, and for its dissipative part. Here we choose the 
simplest L = 0§ — w 2 V 2 L , which models a membrane 
(characterized by its sound velocity v) that can oscil- 
late as a whole (the homogeneous mode) at a frequency 
f^o- This situation could be implemented experimentally 
by attaching the membrane to a compressible substrate. 
Certainly, if the mirror is mounted on a fixed frame the 
displacement field Q (r, t) cannot posses a homogeneous 
mode, but our choice would still be a reasonable approx- 
imation valid for large mirrors as compared to the trans- 
verse size of the optical beam, and low light intensities 
(which induce small deformations in the mirror), simi- 
larly to what happens in cavities with nonlinear elements 
[25|]. With this choice the normal modes are 2D plane 
waves exp(iAc-r), specified by a (transverse) wave vec- 
tor k, with associated frequencies Cl 2 . = CIq + v 2 k 2 , being 
K = \k\. As for the dissipative dynamics we assume the 
simplest case of equal losses 7 K = 7 m for all modes, so 
that using Eqs. (JI|)-([7]) one readily obtains 



d 2 Q + lm d t Q + ( n 2 



v 2 V 2 



Q 



2hk c 



\A\< 



(8) 
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FIG. 2: Illustration of the system's instabilities. We plot the intensity I of the homogeneous steady-state solution as a function 
of the injection E (both quantities are in arbitrary units), with the stable (unstable) regions denoted by solid (dashed) lines. 
We mark with squares both turning points, as well as the Hopf instability leading to oscillatory solutions and the Pitchfork 
instability leading to pattern formation. As explained in the text, three relevant scenarios can be found as the detuning is 
decreased from A = —2 (assuming p > 1 and small mechanical losses): absence of the Pitchfork bifurcation for A above some 
threshold Ao (a), and presence of the Pitchfork bifurcation before (b) or after (c) the Hopf instability. For small rigidities 
(p < 1)) t ne scenario is like (c), but with the Pitchfork bifurcation going to infinity, irrespective of the value of the detuning. 



Equations (|3|) and ((8j) constitute our optomechanical 
cavity model. Before continuing however, and in order 
to give a cleaner presentation, we introduce normalized 
and dimensionless time r = ^ c t and spatial coordinates 
r' = r/l c , together with the variables and parameters 



overbar denotes the homogeneous steady state), with 



Z=^Q, F=tJ 2 ^A, k = Z c K, (9) 



2 2hk c k L tt 7m 

E = jr\ 7 rp £ > = — , 7= — • 
Q V t c oT 7 C 7 C 

Using Eqs. © and © we get 

d T F= (-1 + iA + iV 2 +iZ)F + £, (10a) 
d 2 Z + jd T Z + fig (1 - p 2 V 2 ) Z = fl 2 Q \F\ 2 , (10b) 

where V 2 = I 2 (d 2 + d 2 ) = (d 2 , + <9 2 ,), and we defined 
the 'effective rigidity parameter' 



p = v/Cl l c 



(11) 



so that the larger p the more rigid the mirror behaves. 
Alternatively p is a measure of how local the response of 
the vibrating membrane is to a stress: when p — ¥ the 
response is local while for p — > oo it is completely inte- 
grated, that is, only the homogeneous mode is excited. 
Note finally that in the (just mathematical) limit 7—5-00 
and p — > 0, the mirror field can be adiabatically elimi- 
nated as Z — \F\ 2 and our model reduces to that of the 
self- focusing Kerr cavity 

Homogeneous solutions and their stability. 
When the injected field is a monochromatic plane wave 
propagating along the cavity axis, the amplitude E is a 
constant, which we take as real and positive without loss 
of generality (this just sets the reference phase). In such 
a case (jTUJ) admits spatially homogeneous, steady solu- 
tions, which verify Z = \F\ 2 4= I (in the following an 



E 2 = 



i + (A + /y 



(12) 



The above is the usual state equation of the single- 
mode pendular cavity (or of a Kerr cavity), which pre- 
dicts single-valued (three-valued) response for A > 
(A < — v3)- In the latter case the characteristic I — E 2 
curve exhibits the well known S-shape (Fig. [2]), with the 
turning points located at I± = (-2A ± y/A 2 - 3) /3. In 
the following we will use / instead of E as a control pa- 
rameter because / fixes the steady state and determines 
univocally the value of E through (fT2")) . not the other way 
around. 

This spatially homogeneous steady state can suffer dif- 
ferent instabilities, or bifurcations, leading to the appear- 
ance of new type of solutions, such as time-periodic solu- 
tions or spatial patterns. In the supplemental material we 
explain how to find these instabilities through a standard 
linear stability analysis in which one studies the evolu- 
tion of small perturbations around the steady solution. 
Indeed, we have classified all the possible instabilities by 
using the Routh-Hurwitz criterion ptjj and have found 
that spatial patterns and optomechanical cavity solitons 
can appear under a large variety of situations. However, 
since our goal in this Letter is showing that pattern for- 
mation is possible by means of optomechanical coupling, 
we focus on a region of the parameter space where the in- 
stabilities can be summarized in a simple fashion, leaving 
the exhaustive analysis for a future publication. Hence 
we concentrate in the following on the experimentally 
relevant limit 7 fio (the mechanical oscillations are of 
high-Q) and 7 <g 1 (the damping rate of the mirror os- 
cillations is much smaller than the cavity damping rate). 
As well we restrict our presentation to negative values of 
A (red detuned injection), in particular to A < —2, as 
this is the detuning domain where the analysis is more 
transparent. 
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FIG. 3: (a) Same as Fig. [2]for the parameters of the numerical simulation A = —2.2, 7 = 0.1, fio = 10, and p = 1.13. We have 
found two types of patterns: periodic (hexagonal in 2D) patterns (b), whose average intensity (half the sum of their maximum 
and minimum) is represented as a yellow, dashed-dotted line in their domain of existence; and solitons (c), represented as a 
green, dotted line. Above the Hopf instability we find homogeneous, time-periodic solutions (not represented in the figure) . 



In the considered detuning domain A < — 2 the S- 
shaped response is present. We find that the lower branch 
of the characteristic I — E 2 curve is stable in all cases, 
while the middle branch is intrinsically unstable, as usual. 
As for the upper branch, it is always affected by a ho- 
mogeneous Hopf bifurcation, which occurs at an intra- 
cavity intensity 7hb, the steady state being unstable for 
I > Ihb, see Fig. [2] The expression for Ihb is lengthy, 
but in the considered limit 7 <C 1 we find 



2/hb 



A2+7(l + fi ) 2 /fi2 



(13) 



On the other hand, and this is the interesting point, 
a static, pattern forming instability exists that affects 
the upper branch. It occurs at an intracavity intensity 
Jpb (whose analytic expression is too lengthy to be given 
here) and the steady state is unstable for I + < I < Jpb; 
hence this instability acts oppositely to the Hopf, see Fig. 
[2^b). For p > 1 the instability exists for —00 < A < Ao, 
where Aq is the only real (and negative) solution to 



P 2 = 3- 



2 ^A 2 , - 3 - Ao 



3 + A A + ^Ag-3 



(14) 



Just at A = A the bifurcation is located at the upper 
turning point (Jpe = 1+) and as we depart more from 
resonance (i.e. by making A more negative) it climbs 
through the upper branch, eventually entering the region 
which is already Hopf unstable (Jpb > /hb), see Fig. 
HJc). As the rigidity p is decreased (p — > 1 + ), Jpb — s- 
00 and for p < 1 the whole upper branch is unstable 
against spatial perturbations irrespective of A. Hence 
three basic scenarios exist in the considered limit (7 <C 
f2o, 1; A < —2), depending on the relative position of the 
two bifurcations, as depicted in Fig. [2j 

We conclude that pattern formation is favored by small 
rigidities, as expected, or by large (negative) detunings. 

Numerical integration of the model equations. 
The linear stability analysis so far presented cannot in- 
form us about the type of patterns appearing once the 



bifurcation is crossed, hence a numerical analysis is nec- 
essary. We shall concentrate on the static patterns ap- 
pearing at low intensities (below the Hopf bifurcation), 
under the conditions specified above, leaving the anal- 
ysis of the dynamic patterns for future work. We have 
numerically integrated the model Eqs. (fT0|) in ID and 2D 
by using the split-step method for the special parameter 
choice 7 = 0.1, ft = 10, A = -2.2, and p = 1.13. In this 
case A < Ao = —1.75847, and then the upper branch of 
the steady homogeneous solution has a finite domain of 
stability for I G [Ipb, Jhb] , as in Fig. [2jb). In Fig. [3ja) 
we represent the characteristic curve I — E 2 , together 
with the average intensity corresponding to the different 
patterned solutions we have found. The patterns are of 
two types: periodic patterns, which are hexagons in the 
2D case, see Fig. [3]^b) , and localized structures, see Fig. 
[SJc). We have checked that these localized structures 
are actual cavity solitons (CS): (i) they can be written 
and erased (by using appropriate injected pulses) at any 
point, (ii) when several solitons are present, the deletion 
of one does not affect the rest, and (iii) they can be moved 
across the transverse plane by using, e.g. phase gradients 
of the injected field. 

We would like to remark that in the Ipb > /hb case 
corresponding to Fig. [21c), the pattern formation sce- 
nario is similar, except for the fact that the patterns may 
start pulsing in vicinity of the Hopf bifurcation. 

Conclusions. We have proposed a model for pattern 
formation in an optomechanical system. Our model con- 
siders a cavity formed by a fixed rigid mirror and an oscil- 
lating deformable mirror. We have simplified the descrip- 
tion of the latter by assuming the existence of a homoge- 
neous oscillating mode together with mechanical waves 
(modeled as those of a vibrating membrane) , proving the 
existence of periodic patterns and localized structures. 
Future venues will include the study of the quantum spa- 
tial correlations present on (and between) the optical and 
mechanical fields, which may lead to noncritical squeez- 
ing and multipartite entanglement similarly to what has 
been found in optical parametric oscillators pol . 27 1. 
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Supplemental material 

Derivation of the light field equation. In the fol- 
lowing we deduce equation (|3]) of the main text. To this 
aim we use the approach of reference |22j , which consists 
in propagating the complex amplitudes A± (z, r, t) along 
a full cavity roundtrip, for what we assume that they is 
slowly varying in space and time so that they satisfies 
the paraxial wave equation 



(d z ± c~ l d t ) A ± 



i ■ 



(15) 



Given the amplitude A+ (z = L,r,t), after reflection 
on the deformable mirror we get A- (L, r, t) e~ lfeLL = 
- exp [2ifc L Q (r, t)} A + (L, r, t) e ifcLi where Q (r, t) repre- 
sents the displacement of the mirror from its rest position 
(Q = at rest). This amplitude A_ [L, r, t) propagates 
from z = L to z = giving rise to a new amplitude 
A_ (0,r,t+ |i c ) = U L A_ (L,r,t), where t c = 2L/c is 
the cavity roundtrip time and Ul = exp [i(L/2fcL)V 2 ] 
corresponds to the paraxial propagation operator in free 
space. After reflection onto the coupling mirror a new 
amplitude A + (0,r,t+ \t c ) = -y/RA- (0,r,t+ \t c ) + 
V^Ajnj (0,r,t+ \t c ) is got, with R and T the reflec- 
tivity and transmissivity factors of the coupling mir- 
ror, respectively (i? + T = 1 is assumed: lossless mir- 
ror). Finally propagation from z = to z = L yields 
A + (L, r, t + t c ) = U L A + (0, r, t + \t c ) . Adding all parts 
together one gets 



A(r,t + t c ) = VRe 2ikLL Ul exp [2ifc L Q (r, t)} A (r, t) 



VfA ini (L,r,t + t c ), 



(16) 



where we used U^A^ (0, r,t+ ht c ) = A in j [L, r,t + t c ). 
We now take into account that R —> 1 so that \/H = 
y/l — T — >■ 1 — ifT: Remember that T is very small. 
Next we assume that light is almost resonant with the 
cavity, specifically we impose that 2 (cjl — w c ) L/c = 5 is 
of order T, where w c is the cavity longitudinal mode fre- 
quency (hence uj c = mirc/L, m £ N) closest to wl, what 
allows approximating exp(2ifcL^) = cxp(2iiOhL/c) ~ 
1 + iS. We assume as well that k^Q (r,t) is of order T 
(the mirror displacement/deformations are much smaller 
that the optical wavelength), so that exp [2ik^Q (r,t)] w 



1 + 2iki,Q(r, t). Similarly we assume that the effect 
of diffraction is small (this implies that both mirrors 
must be sufficiently close each other, either physically 
or by means of lenses) so that we can expand XJ\ ~ 
1 + \(Ljkx)V\. All these assumptions imply that the 
overall variation of A between consecutive roundtrips 
is very small and then one can approximate df A by 
[A (r, t + t c ) — A (r, t)] i" 1 . With all these approxima- 
tions we get, to the lowest nontrivial order, 



dtA (r, t) = 7c I -1 + iA + U 2 C Vl + i 



AkL 
T 



A + 7c £, 
(17) 

where 7 C = T/2t c = cT/AL is the cavity damping rate, 
A = (wl — uj c ) /7 C is the dimcnsionless detuning parame- 
ter, I 2 — 2i/fc L T is the diffraction length (squared), and 
£ (r, t) = 2T~ 1 / 2 A- m - i (L, r, t + t c ) is a scaled version of 
the injection field amplitude. 

Details of the linear stability analysis. In the 
following we explain how we have performed the sta- 
bility analysis of the homogeneous, stationary solution 
(|T2"T) given in the main text. We have followed the stan- 
dard linear stability analysis, by studying the evolution 
of small perturbations (SF, 5Z) added to the steady solu- 
tion (F, Z) . Upon linearizing the model equations (p~0|) 
with respect to (SF,SZ), expressing them in terms of 
the normal mode basis as <5F(r',r) = ^ k 0k (t) e lk ' r 
and 5Z{v',t) = £ k Ck (r) e ik r ' (note that £ k = Ck 
because SZ is a real field), and equating coefficients 
of like exponentials, we get, for each k, a linear sys- 
tem of differential equations for the modal perturbations 
v = (</>ki Ck) ■ Owed to this linearity and the time 
invariance of the system, its solutions have the form 
v ( T ) — v (0) &Xr ■ Upon making such substitution a ho- 
mogeneous linear system of algebraic equations in v (0) 
is got, whose condition for existence of nontrivial solu- 
tions reads C (k 2 ; A) = J2t=o Cn (^ 2 ) ^™ = 0> where k is 
the modulus of k. After simple algebra we get ca = 1, 
c 3 = 2 + 7, and 



C2 
Cl 

CO 



i + 2 7 + n 2 + (A fc + If 



1 



21 A k 



(A fe + /y 



1 



(18a) 
(18b) 

(18c) 



where fi? = 



n 2 



1 + p 2 fc 2 ) and A fe = A — k 2 . We observe 
that the growth exponents A, solutions to C (fc 2 ; A) = 0, 
depend on k and not on k because of the rotational invari- 
ance of both the steady state and the model equations. 
Whenever Re{A} < for all k the steady state is sta- 
ble, while if Re{A} > for some k it is unstable. The 
condition Re{A} = thus defines a possible instability, 
or bifurcation, which is met either when A = (static 
instability: Co = 0) or when A = \\fc\Jcz (self-pulsing, 
or Hopf instability: C1C2C3 = C4C 2 + c 2 co). On the other 
hand, when the bifurcation is associated with k = the 
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new state is spatially uniform (homogeneous instability), 
while if k ^ the instability is pattern forming. 
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